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Abstract 

In this paper, we present a theoretical justification of the 0-1 test for chaos. 
In particular, we show that with probability one, the test yields for periodic 
and quasiperiodic dynamics, and 1 for sufficiently chaotic dynamics. 

1 Introduction 

In [5], we introduced a new method of detecting chaos in deterministic dynamical 
system in the form of a binary test. The method applies directly to the time series 
data and does not require phase space reconstruction. As explained in [5], with 
probability one the test gives the output = for quasiperiodic dynamics and 
K = 1 for sufficiently chaotic dynamics. 

In [6], we proposed a simplified version of the test that is more effective for sys- 
tems with a moderate amount of noise. The effectiveness of the new method was 
demonstrated for higher- dimensional systems in [6] and for experimental data [3]. 

The main aim of this paper is to put the simplified version of the test on a rigorous 
footing, going far beyond the results indicated in [5] for the original test. In addition, 
our analysis of the test leads to a significant improvement which was used in our 
paper [7] detaihng the implementation of the test. 

We first recall the simplified form of the test proposed in [6]. Let f : X —>■ X 
be a map with invariant ergodic probability measure fi. Let f : A M be a scalar 
square-integrable observable. Choose c G (0,27r), x E X, and define 

n-l 

Pe(n)=5^e^^'=t;(fx). (1.1) 
i=o 
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Next, define the mean-square displacement 

1 ^ 

Me(n) = lim - V M + n)- p,{j)\\ (1.2) 

Finally, let 

= li„, ^MM. (1.3) 

N^oo log n 

The claim in [6j, substantiated in this paper, is that typically (i) the limit Kc exists, 
(ii) Kc G {0, 1}, and (iii) Kc = signifies regular dynamics while Kc = I signifies 
chaotic dynamics. 

Remark 1.1 (a) The definition of Pc{n) in (11.11) is slightly different from in [3l O [7] 
where pdn) = ^JZq cos(jc)f (/%). In the current paper it is natural to simplify 
analytic calculations rather than numerical computations, but apart from that the 
methods are equivalent. 

(b) For fixed c, it follows from the ergodic theorem that the limit Mc{n) in (11.21) 
exists for almost every initial condition x and the limit is independent of x. The 
common limit is 

Mcin)= I \pc{n)\''d^l=\\pc{n)\\l 
Jx 

To see this, compute that pc{i + n) — Pc{j) = e^^^pdn) o and so Mc{n) = 
limTv^oo XljL'o^ bc(^)P o P which converges to the space average Jj^ \pc{n)\'^dfi al- 
most everywhere. 

(c) Strictly speaking, the limit Kc in (11.31) need not be well-defined. Of course, 
K^ = limsup^^QQ log Mc(n)/ log n is well-defined, and it follows from Proposition II. 41 
that K^ G [0, 2] for all c. (In the case of periodic dynamics, Kc = 2 for isolated values 
of c.) 

Example 1.2 Consider the logistic map / : [0, 1] — > [0, 1] given by f{x) = ax{l — x) 
for < a < 4. This family of maps is particularly well-understood [121 1]: we can 
decompose the parameter interval according to [0,4] = V U C U J\f where J\f has 
Lebesgue measure zero and the asymptotic dynamics consists of a periodic attractor 
(of period g > 1) for a G P and a strongly chaotic attractor consisting of g > 1 
disjoint intervals for a G C (satisfying the Collet- Eckmann condition). 
We obtain the following result: 

Proposition 1.3 Let v : [0,1] ^ W be Holder. Let a E V UC and define q as above, 
(a) If a eV , then Kc = for all c ^ li^j jq. 



2 



(b) If a E C, then = 1 for all c 7^ ^nj/q unless v is infinitely degeneratqj. 
Hence, the test succeeds with probabihty one for logistic map dynamics. 

Part (a) holds for general periodic dynamics (and all continuous observables) . In 
Section [21 we prove that the test yields Kc = 0, for almost all c, for quasiperiodic 
dynamics, provided we make smoothness assumptions on v. This justifies our claim 
that Kc = for regular dynamics. 

The chaotic case is discussed extensively in Section |3l In particular we obtain 
Kc = 1 under various assumptions: 

(i) Positivity of power spectra 

(ii) Exponential decay of autocorrelations 

(ill) Summable decay of autocorrelations plus hyperbolicity 

(In fact, (ii) and (iii) are sufficient conditions for (i).) 

In many situations, including the logistic map with a G C, it is necessary to 
consider instead of /, and autocorrelations decay only up to a finite cycle (of 
length q). As shown in Section [3l criteria (ii) and (iii) generalise to this situation. 



Summable decay without hyperbolicity assumptions Without making hy- 
perbolicity assumptions, we have no definitive results when autocorrelations decay 
subexponentially. However, there is some partial information discussed in Section |H 
If the autocorrelation function is summable, then the power spectrum S{c) exists for 
all c G (0,27r) by the Wiener-Khintchine Theorem [9j, implying that 

Mcin) = S{c)n + o{n). 

Under slightly stronger assumptions on the decay rate 

Mc{n) = S{c)n + 0{1). 

In the former case, K^ = limsup^^^^^ logMc(n)/logn G [0, 1]. In the latter case, Kc 
exists and takes the value or 1 depending on where S{c) = or S{c) > (but see 
Remark 13.71) . 

Again, we obtain similar results if autocorrelations are summable up to a finite 
cycle. 

^ Lying in a closed subspace of infinite codimension in the space of Holder functions 
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Improved diagnostic in the test for chaos The o{n) and 0(1) terms above are 
nonuniform in c but in Section S] we show that the source of nonuniformity is easily 
dealt with. Define 

D,{n) = M,{n) - {Evf- . 

1 — cos c 

Here Ev = f^vdfi denotes expectation with respect to fi. Under the above conditions 
we obtain Ddn) = S{c)n + o{n) (hence = limsup^^^^ log Mc(n)/ log n G [0,1]) 
for summable autocorrelation functions and Dc{n) = S{c)n + 0(1) (hence Kc takes 
the values either or 1) under slightly stronger conditions on the decay of the au- 
tocorrelation function as before, but the o{n) and 0(1) terms are now uniform in c 
(see Section HI). In |7j, we proposed using Ddn) instead of Mdn) in the numerical 
implementation of the 0-1 test, and demonstrated the improved performance of the 
test. 

Nonsummable decay The summability condition in the Wiener-Khintchine Theo- 
rem can be weakened considerably. For example, if autocorrelations decay at a square 
summable rate (including k~'^ for any d> then the power spectrum exists almost 
everywhere and so Mc(n) = S{c)n + o{n) for almost every c. (In this generality there 
is no uniformity in the error term for Dc{n.) This and related results is discussed in 
Section [51 

Correlation method Our emphasis in this paper is on understanding the proper- 
ties of the limit as defined in (11.31) . However, in [7], we proposed computing 
as the correlation of the mean-square displacement Mc{n) (or Dc{n)) with n. The 
advantages of this approach were demonstrated in [7j. In Section [HI we verify that the 
theoretical value of Kc remains for regular dynamics and 1 for chaotic dynamics. 

The paper concludes with a discussion section (Section [7]). We end the introduction 
by proving that G [0, 2] as claimed in Remark 11.1( c). 

Proposition 1.4 Let = limsup^^o^ log Mc(r;,)/ log ra. Ifv is not identically zero, 
then K+{c) G [0,2] for all c. 

Proof By definition, ||pc(^)||2 < ''^ll'^lh so that < Mcin) < n^\\v\\l. Hence < 2. 

To prove the lower bound, we use the fact that [It'll 2 > 0. It suffices to show that 
limsup„^oo Mc{n) > for each fixed c. Observe that pcip- + 1) = e^'^f o /" -|- pc(^) 
so that 

\\Pc{n + 1)||2 > We'^^'v o /"i|2 - \\pc{n)\\2 = \\vh - \\Pc{n)h. 

Hence < ||f II2 < ||Pc(^)||2 + ||Pc(^ + 1)||2- It follows that ||pc(^)||2 7^ 0, and so 
lim sup„^oo Mc(n) > as required. I 
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2 The case of regular dynamics 

Part (a) of Proposition II .31 is a simple direct calculation. If / : X ^ X is a map with 
a periodic orbit of period q and f : X — M is continuous, then we obtain Kc = ^ for 
all c 7^ 271] /q. (For isolated resonant values c = l-nijq a simple argument using the 
Fourier series for v shows that typically "Pcip) will grow linearly implying i^c = 2.) 
In the case of quasiperiodic dynamics, we require additional smoothness assumptions 
on the observable v. The test then succeeds with probability one. 

Theorem 2.1 Suppose that X = T™ = M'"/(27rZ)™ and that f : X X is given by 
f{x) = X + uo mod 27r. //f : X M zs C" with r > m, then Kc = for almost every 
c e [0,27r]. 



Proof Write f : X ^ M as a m-dimensional Fourier series v{x) = ^^g^m f 
where = ve- Then 



i£-x 



n— 1 n— 1 



^ j=0 p,q=0 "^^ 

n— 1 „ 

J2 {n-\k\)e^^ j vo f\^\vdx = si + --- + sn, (2.1) 



k=-{n-l) 

where 



m— 1 „ 



j=-(m-l) 

We show that Mdji) is bounded (as a function of n) for almost all c. Compute 
(formally) that 



n-l 



i=-(n-l) ^ 

|ti^|^(e*^'^^^''^'' — l)"^(e*"'^'^"'"^''^^ — e"**^""-^^'^'^"'"^''^^). (2.2) 



Hence 



E, ,9 1 — cosn(c + I ■ Oj) , , 

^ 1 — cos(c + ^ ■ C<j) 

It remains to show that the series (12.21) . (12. 3p converge. We may ignore the £ = 
term in these series (these terms are obviously bounded in n). The smoothness 
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assumption on v implies that \ve\ = 0{\i\ Let e > 0. For almost every c G (0, 2tt) 
there is a constant do > such that 

dist \c + e-uj, 2itZ\ > do\i\-^"'+'\ (2.4) 
for all £ G Z'" - {0} (cf. fTB]). Hence |e*("+^-^) - 1| > ^11^1"^™+'^ and so 

oo oo 

e&Z"^-{o} k=i \e\=k k=i 

oo 

= C5^A;-(W-m-e/2)) 

k=l 

provided we choose e > so small that r > m + e/2. This shows that (12. 2p converges 
and is bounded independent of n, and similarly for (12.31) . I 

Remark 2.2 The extra smoothness of v is required to circumvent the small divisor 
problems associated with quasiperiodic dynamics. We also require a Diophantine 
condition on c, satisfied by almost every c G [0,27r]. However, there is no restriction 
on uj. 



3 The case of chaotic dynamics 

It is our intention to show that = 1 for all almost all c (and reasonable observables 
v) for sufficiently chaotic dynamical systems. We proceed along three distinct but 
related avenues, all of which extend Example 11.21 of the logistic map: (i) positivity 
of power spectra; (ii) decay of autocorrelation functions; (iii) hyperbolicity of the 
dynamical system. 

Recall that for a square-integrable observable v : X ^ W the power spectrum 
S : [0, 27r] [0, oo) is defined (assuming it exists) to be the square of the Fourier 
amplitudes of f o per unit tim^, and is given by 

, „ 71-1 ^ 

S{c) = lim - / o f^^ dfx = lim -Mc(n). 

n^oo n Iy ^ — ^ n— >oo n 

In other words, Mc{n) = S{c)n + o{n). The following result is immediate: 

Proposition 3.1 Let c G [0,27r]. Suppose that S{c) is well-defined and strictly posi- 
tive. Then Kc = 1. I 

^Often e'^" is replaced by ^'^i-i^h'^ in the literature, but this is just a rescaling of the domain. 
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In particular, if the power spectrum is well-defined and positive almost everywhere, 
then we obtain = 1 with probability one. 

Next, we consider the autocorrelation function p : N ^ M. given by 

p{k) = j vof^vdp—{^j vdp^ . 

This is well-defined for all observables v. 

If p{k) is summable (i.e. YlT=o\Pi^)\ ^ then it follows from the Wiener- 
Khintchine theorem [9] that for c G (0, 27r), 

oo 

Sic) = e^'^p{\k\). 

k=—oo 

Note that the right-hand-side defines a continuous function on [0,27r]. 

Proposition 3.2 Suppose that v : X ^ M. lies in L'^{X) and v is not constant 
(almost everywhere). If the autocorrelation function p{k) decays exponentiall^, then 
Kc = 1 except for at most finitely many choices of c & [0, 2tt] . 

Proof Since p{k) decays exponentially, g{c) = YlT=-oo e^'^Vd^l) is analytic on [0, 2n]. 
Since v is not constant, p(0) = f^v'^dp — {fj^vdpY = f^{v — J^v^dp > 0, and 
hence g is not the zero function. By analyticity, S{c) = g{c) > except for at most 
finitely many values of c and hence Kc = 1 except for these values of c. ■ 



Decay of autocorrelations up to a finite cycle Recall that / : X ^ X is mixing 
if lim^^oo p{k) — >■ for every observable f : X — > M. The system is mixing up to a 
finite cycle (of length g > 1) if X = Xi U ■ ■ ■ U X^ where /(X^) C X^+i (computing 
indices modg) and : Xj Xj is mixing (with respect to Pg = qp\Xj) for each 

If g > 2, then decay of autocorrelations holds only for degenerate observables. 
The natural property to require is exponential decay for f^. Given an observable 
f : X — * M, define for j = 1, . . . , g and m = 0, . . . , g — 1, 

pvof^,v,Xj [kg) = V o f"" o f''" V dpj - vo dpj / v dpj. 

J Xj J Xj J Xj 

Definition 3.3 The autocorrelations of v are summable up to a q cycle if for each 
j = 1, . . . , g and m = 0, . . . , g — 1, the series J2T=o \Pvof"^,v,Xj{kq)\ is convergent. 

The autocorrelations of v decay exponentially up to a q cycle if Pvof"^,v,Xj{kq) 
decays exponentially as A; — > cxd for each j = 1, . . . , g and m = 0, . . . , g — 1. 

^There exist constants C > 1, r e (0, 1) such that \p{k)\ < Ct'' 
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Theorem 3.4 If the autocorrelations of v are summable up to a q cycle, then 

oo 

S{c)= e'^'^gr, for all c ^ 2-nj/q, 

r=—OD 

where writing r = kg + m with A; G Z and m G {0, 1, . . . , g — 1}, 

gr = gkq+m = Pvof^,v,Xj{kq). 

i=i 

Proof Define 

9-1 



e=o ^ 1=1 "^-^j 



Vq,c O f'"^Vg^c dfXj 



Vq,c 



By [lii Theorem A. 2], 



S{c) = J2 ^'"''P.Ak), 



k=—oo 



for c 7^ 271] /q. Compute that 



1 /" 

1 ^"^X 



i=i -^^^ 



Vq,c dpj 



e'^'-'^\j vof^of'^vof^'df^,- I vof^df^.l vof^'d^^, 

q-1 q 

fkq 



V dfij+e> - / vo f d^j. 



1 V % ^ 

Q ^ ^ 

^ j = l £/'=( 

^ ^/'=0 j = l 

^ ^/'=0 j=l "^-^J 

-E E vof^of'^'^vdfi,- vof^dfij vdfi,>j 



3+1' 

t> o f^~^ dfij / f c//ij 



^3 I " 

IX-i 



3=1 s=-{q-l) 
q q-l 



^ i=l 5=-(g-l) 



For r = kq + m we obtain at most two nonzero contributions to gr, namely s = m 
in Pq^c{k) and s = —{q — m) in Pq^c{k + !)• Hence 



gr = - 



(kq + q))=J2 

(kq) 



as required. 



Corollary 3.5 If the autocorrelations of v are summable up to a q cycle, then S{c) 
exists and is continuous except for removable singularities at c = 2'irj/q. 

If the autocorrelations ofv decay exponentially up to a q cycle, then S{c) is analytic 
except for removable singularities at c = 27c j/q. If moreover v\Xj is not constant 
(almost everywhere) for at least one j, then Kc = 1 except for at most finitely many 
values of a. 

Proof The statements about continuity and analyticity are immediate from Theo- 
rem I3.4[ In particular, if there is exponential decay up to a g cycle, then the function 
g{c) = Yl'^-oo ^^^"^dr is analytic and hence nonzero except at finitely many points 
provided go 7^ 0. If on the other hand, 5^0 = 0, then 



Remark 3.6 Definition 13.31 and Theorem 13.41 are significant improvements on the 
corresponding material in [T^ Appendix]. 

Remark 3.7 We have seen that exponential decay of autocorrelations (up to a g 
cycle) guarantees that Kc = 1 with probability one. Surprisingly, it seems nontrivial 
to weaken the exponential decay hypothesis. The proof of Proposition 13.21 relies 
crucially on analyticity of the power spectrum. Even if we assume sufficiently rapid 
decay that g{c) = XlfcL-oo ^^'^Vd^D is C°°, then we face the difficulty that the only 
restriction on the zero set of a function is that it is a closed set. 

Suppose that summable decay of correlations holds for a large class of observables 
B with the property that there is an interval / C (0, 2it) such that the Fourier series 
g{c) is identically zero on / for all v ^ B. The proof of Proposition 13.21 shows that 
for every nonconstant observable v & B, there is an interval J C (0, 2-7?) on which 
g{c) > on J. For such examples, where the power spectrum vanishes on an interval 
/ and is typically positive on an interval J, the 0-1 test is inconclusive: we obtain 



= ^?o 




so var(f |Xj) 



0, and hence v\Xj is constant, for each j. 
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Kc = 1 with positive probability for nonconstant observables in B, but for all f G i3 
there is a positive probability that G [0, 1) or that Kc does not even exist. This 
situation seems highly pathological, but we do not see how to rule this out. 

Hyperbolicity We can overcome the unsatisfactory aspects of Remark 13.71 by as- 
suming some hyperbolicity. In the Collet-Eckmann case (a G C) for the logistic map, 
it is known that Holder observables enjoy exponential decay of correlations up to a 
finite cycle, so we can apply Theorem 13.41 Alternatively, [H] shows that the power 
spectrum is bounded away from zero for all c 7^ 27ij/q, so we can apply Proposi- 
tion 13.11 These comments apply to all maps in the following classes: 

• Uniformly expanding maps; Uniformly hyperbolic (Axiom A) diffeomorphisms. 

• Nonuniformly expanding/hyperbolic systems in the sense of Young [17], mod- 
elled by a Young tower with exponential tails. These enjoy exponential decay 
of correlations (up to a finite cycle) for Holder observables. This covers large 
classes of dynamical systems, including Henon-like maps, logistic maps and 
more generally multimodal maps satisfying Collet-Eckmann conditions, and 
one-dimensional maps with Lorenz-like singularities [2]. 

Young [18] weakens the decay rates assumed for tower models for nonuniformly 
expanding/hyperbolic systems. Holder observables now have subexponential decay 
of correlations (up to a finite cycle). Provided the decay rate is summable, the 
argument of [H] still applies, and S{c) is bounded away from zero (except for infinitely 
degenerate observables) . 

Example 3.8 A prototypical family of examples is the Pomeau-Manneville inter- 

{x(l + 2°x"); < a; < - 
~ ~ ^ where 
2x — 1; ^ < X < 1 

a G [0, 1) is a parameter [T6l [TT] . When a = this is the doubling map with expo- 
nential decay of correlations for Holder observables, so Proposition 13.21 applies. For 
a > 0, let /? = ^ — 1. Then decay of correlations for Holder observables is at the rate 
[8]. By [H], the power spectrum is bounded below in the summable case (/3 > 1, 
equivalently a < |) and so = 1 for all c G (0, 2tt). 

Remark 3.9 Numerical experiments for intermittency maps indicate that (i) Kc = 1 
and (ii) the power spectrum exists and is bounded below, even in the nonsummable 
case P < I, equivalently a G [|, 1). It remains an interesting problem to prove these 
statements. By Corollary 15.3 j we are at least assured that the power spectrum exists 
(and hence Kc = 1 with positive probability) for /5 > | (a < |). 
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4 Summable decay without hyperbolicity 



Recall that if the autocorrelation function is summable (i.e. YlkLi\pi^)\ < 
then the power spectrum S{c) exists and is continuous for all c e (0, 27r). Indeed, 
'S'(c) = X]fcL-oo ^^'^Vd^D oil (0,27r) by the Wiener-Khintchine Theorem P). Hence 
Mc{n) = S{c)n + o{n). As noted in Remark 13. 7^ S{c) is not identically zero for 
nonconstant observables, and hence there is an interval of values of c for which S{c) > 
0. In particular, Kc = 1 with positive probability. 

In this section, we discuss the error term o{n) in more detail. As mentioned in 
the introduction, this leads to the improved diagnostic Dc{n) for chaos used in |7]. 

We begin with a formal calculation to express the mean square displacement as 
follows: 

Proposition 4.1 M^in) = EL-n(^ - \k\)e''"'p{\k\) + {Evf 
Proof First, note that 

n-l 



Me(n) = I \Pc{n)\^ d^^=Y^ (^'^"'"^^ [ v o fP v o f dfi 

•^^ p,q=0 

n-l „ 

= e'^^""^' / V o f\P-'i\ V d^l 

p,q=0 ^ 
n-l 

p,g=0 

n n-l n-l 

= 5^ (n - |A:|)e^'=^p(|A:|) + {EvfY^^^P^H^ 



p=0 q=0 



Finally 



n-l n—1 



^ _ g«nc _ ^ mc _ 



e > e 



7 o — : : — . 

/L^ I _ gtc I _ g-?.c I _ g 

p=0 g=0 g 

The second term in the expression for Mc{n) is bounded in n for fixed c, but is 
nonuniform in c. Since the term is explicit, it is convenient to remove it. (As demon- 
strated in [7], this is also greatly advantageous for the numerical implementation of 
the test.) Hence we define 

D^{n) = M^{n) - (Evfil - cosnc)/(l - cose) . 
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Then Mc{n) = Dc{n) + 0(1), so it suffices to work with Dc{n) from now on. By 
Proposition I4.H 

n 

D^{n)= J^i^-lkDe^'^pm). (4.1) 

k=—n 

Theorem 4.2 Suppose that p{k) is summahle (i.e. YlT=i lp(^)l <^ Then for all 
ce (0,27r), 

Dciji) = S{c)n + e(c, n), 

where 

oo n 

\e{c,n)\<2n |p(A;)| + 2 ^ A;|p(A:)| = o(n). 

fc=n+l fc=l 

/n particular Dc{n) = S{c)n + 0(72) uniformly in c. 
Proof Write 

n 

D,{n) =J2{n- \k\)e'^'p{\k\) = S{c)n + e(c, n), 

k=—n 

where 

00 n 

e(c, n) = -n (^'^^ + e"''"')p(A;) - J2 ^^^'^^ + e"''"')p(A;). 

fc=n+l k=X 

It remains to show that X]fc=i ^|p(^)l = o{n). Let s„ = YMt=i |p(^)l and let L = hms„. 
Define the Cesaro average (T„ = }[Y^^=i^k, so L = hm(T„. Then ^J2k=i^\pi^)\ ~ 

Under stronger assumptions on the decay rate of the autocorrelation function p{k), 
improved estimates for the o{n) term are available. 

Theorem 4.3 Suppose that Yl'k=i ^lp(^)l < Then 

Dc{n) = S{c)n + 6*0(0) + e(c, n), 

where So{c) = —J2T=-oo^^^^\^\p(\^\) '^^ continuous on [0,27r], S{c) is on [0,27r], 
and 

00 

\e{c,n)\<2 J2 {k-nMk)\=o{l). 

k=n+l 

In particular, D^{n) = S{c)n + 0(1) uniformly in c. (Hence G {0, 1} for all 
ce (0,27r).; 
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Proof Compute that 

n 

D,{n) = 5^ (n - IkDe^^'pilkl) = S{c) n + S^{c) + e(c, n), 

k=—n 

where e{c,n) = YlT=n+ii^ ~ ~^ ^~^^'^)pi^)- ' 

Corollary 4.4 Suppose that \p{k)\ < Ck~'^ for k > 1. If 1 < d < 2, then in 
Theorem \4-^ 



If d>2, then in Theorem\4.3 



\e{c,n)\<2C{-————+^^ ^ 



(rf- 2) n) 

Proof The first term of e(c, n) in Theorem 14.21 is dominated by 

°° /-co 

2Cn k-'^ <2Cn / a;"^ < 2Cn-('^-2V(^^ - !)• 



fc=n+l 

The second term is dominated by 



2Cy^k^"^ <2c(l+ [ x^-Ux) =2Cn~^'^-^^/{2-d)-2C{d-l)/{2-d) 
k=i ^ ' 

< 2Cn-('^-2)/(2-rf). 



Combining these terms gives the result for 1 < d < 2. 
If d > 2, then by Theorem] 



\e[c, n ) 



oo 

d 



noo POO 

<2C x^-'^dx-2Cn( x''^ dx - n 



d-2 d-lJn'^-'^ n'^-^y 
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5 Nonsummable decay of correlations 

In this section, we reformulate the 0-1 test in terms of Cesaro averages, and give 
surprisingly weak sufficient conditions under which S{c) = limn^oo fx\Pc{'>^)\'^df^ 
exists (for typical values of c). 

Let ttfc G C be a sequence with partial sums Sn = Ylk=-n and set cr„ = 
nX]fc=o"^fc- Recall that the sequence is Cesaro summahle if lim„_,ooO'n exists. 
If hm„^oo Sn ^ L then lim„^oo o'n = L (the converse is not true). 

Defining s„ and ct„ as above with = e'"'"^ p{\k\) , we obtain 

1 " 1 
^n = - $^(n-|A:|)p(|A:|)e^^^=-D,(n), 

k=—n 

the last equality following from (14. ip . Since S{c) = lim„^oo ^Dc{n), we have proved 
the following result. 

Lemma 5.1 Let c G (0,27r). Suppose that the sequence au = e^'"^p{k) is Cesaro 
summahle with limit L{c). Then S{c) = L{c). In particular, if L{c) > 0, then 
K, = l. I 

By Fejer's theorem [TOj, a special case is provided when p(|fc|) are Fourier coeffi- 
cients of an integrable function. 

Theorem 5.2 Suppose that p{\k\) are the Fourier coefficients of an L} function 
g : [0,27r] —>■ M.-^. Then the sequence e^^^p{k) is Cesaro summahle to g{c) al- 
most everywhere. In particular, Lemma \5. 1\ holds for almost every c G (0,27r), with 
L{c)=g{c). 

If g is continuous, then the convergence is uniform in c ( and holds for every 
c G (0, 27r) ). In particular, Dc{n) = S{c)n + o{n) uniformly in c. 

Proof We have written ct„ = ^-Dc(n) = ZlL-n(l ~ ^)e*^''p(l^l)- This is an-i{g,c) 
in [ini p.l2 (2.9)]). 

If g is continuous, then by Fejer's Theorem (|TOl Theorem 2.12]), (Jn ^ g uniformly, 
and so Dc{n) = nan = ng{c) + o{n) uniformly in c. 

For general g E L^, it follows from the discussion in (TUl pp. 19-20] that an ^ g 
almost everywhere, so that Dc{n) = nan{g, c) = ng{c) + o(n) for almost every c. ■ 

Corollary 5.3 Suppose that J2k>iPW^ < Then the first statement of Theo- 
rem \5.^ applies with g{c) = YlT=-oo^^^'^Pi\^\) ■ 

In particular, S{c) exists almost everywhere, and Dc{n) = S{c)n + o{n). 
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Proof Since X]fc>iP(^)^ < ^5 the function g{c) = X]fel-oo ^^''Vd^D lies in 
and the Fourier coefficients of g are precisely p(|A:|). Hence, we can apply the ffist 
statement of Theorem I5.2[ ■ 



6 Correlation method 

In [Tj, we proposed computing Kc as the correlation of the mean-square displace- 
ment Mc(n) (or Dc{n)) with n, rather than computing the limit of log M(n)/ log 
as in (11.31) . In this section we verify that the theoretical value of Kc remains for 
regular dynamics and 1 for chaotic dynamics. 
Given vectors x, y of length n, we define 

^ n 1 " 

cov{x, y) = -y^{x{j) - x){y{j) - y), where x = -V'x(j), 
n ^ — ^ n ^ — ^ 

i=i j=i 

var(x) = cov(x, x) . 

Form the vectors ^ = (1, 2, . . . , n) and A = {Dc{l), -Dc(2), . . . , Dciji)). (In particular, 
var(^) = j^(^^ — 1).) Define the correlation coefficient 

Kc = hm corr(e. A) = G [-1, 1] . 

A/var(^)var(A) 

6.1 Quasiperiodic case 

For quasiperiodic dynamics, we have the following analogue of Theorem 12.11 However, 
we require stronger regularity for the observable v, and a Diophantine condition on 
the frequency u (in addition to the condition on c). 

Theorem 6.1 Suppose that X = T™ = M'"/(27rZ)™ and that f : X ^ X is given 
by f{x) = X + u mod 2tt. Let v : X —>■ be a nonvanishing observable with 
r > 3m/2. If Kc is computed using the correlation method, then for almost every 
uj G [0, 2tt] we obtain Kc = for almost every c G [0, 2tt]. 

Proof The proof of Theorem 12.11 shows that Mc{n) (equivalently Dc{n)) is bounded 
for almost every c provided r > m. We show for r > 3m/2 that var(A) = 0(1) 
and cov(^. A) = 0(1) for almost every u and c. Moreover, we use the fact that v is 
nonvanishing to show that var(A) = a + 0(1/?!), where a > 0. It then follows that 

corr(^. A) = ^ ^^ ^ = 0(l/n), 

^in^-l)^/a + 0{l/n) 
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as required. 

The starting point is the calculations fl2.ip and fl2.3l) which gives 



Mc{3) = 1 / ^ , / =J2Mcosjee - 1 , 

^ 1 — cos(c + t ■ u) ^ 

where = c + £ ■ cj, W£ = — |t>£p(l — cos Qs)^^ ■ Hence 

^c(j) = ^ W£(cos j6'^ - 1) - C{\ - cos jc), 

where C = {EvY/{l — cose). Adding a constant (independent of j) to Dc{j) does 
not alter the value of corr(^, A) so we may replace Dc{j) by 

-^c(j) = Y'^i cosjOi + C cos JC, 
e 

when proving that cov(^, A) = 0(1) and var(A) = an + 0(1) with a > 0. Hence it 
suffices to show that for almost every uj and c there exists a > such that 

n n n 

5^5,(j)=0(l), 5^j5,(j) = 0(n), 5^5,(jr = an + 0(l). 

j=i j=i j=i 

Formally, 

n n n. 

cos jc, (6.1) 

j=i i j=i j=i 

n n n 

i=i ^ j=i i=i 

n n 

E -Oc(jT = E E co^i^^ cos j6'f' 

i=i 3=1 

n n 

+ 2C Wi E^ cos jccosjOg + cos^ jc. (6.3) 



For y9 G (0, 2tt), we have 



y e^j'p = e^^- i, y je'^^ = — . 
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In particular, ^"^j^cosjV = 0(1) and X]y=ii'^osjV = 0{n). By fl6.ll) and fl6.2p . 



formally we have ^YTj=\ ^c{j) = 0(1) and J2^=ijDc{j) = 0{n). Turning to (16.31) . 



n 1 " 1 

^cos2jc= -^(l + cos2jc) = -n + 0{l) 



2 

3=1 3=1 



for n and formally 

n 1 

^ W£^cosjccosj6'^ = - ^ (cos j(c + 61^) + cos j(c - Qt)) = 0(1), 

I 3=1 I i=i 

while 

n \ ^ 

^ cos j6'£ COS j6'£/ = 2 XI ^{'^'^^^{(^e + Oe) + cos j{9i - 9e)} 
j=i £/' j=i 

11*^ 1 " 

= - X w^n + 2 X] cos 2j6'^ + 9 X]"lC°^-^*^^^ + 6'f') + cos j(^£ - 6'f')} 

= i^^,^n + 0(l). 



2 

j=i e^fj j=i 



Hence X]j=i ^c{j)^ = an + 0(1), with a = |(^^ w| + O^). If f is nonvanishing, then 
Vi, and hence W£, is nonzero for at least one i so that a > 0. 

It remains to justify the formal calculations. We give the details for the first term 
in (16.31) focusing on the most difficult expression 

n 

I = WjWi' ^{cos j{9e + 9i>) + cos j{9e - 6*^/)}. 
li^i' j=i 

Let e > 0. We assume the Diophantine conditions (12. 4p and 

dist |2c + £-cu,27rZ| > rfal^r^""*"^ dist |£ ■ u;, 27rZ| > c/sl^r^"^'^ (6.4) 



17 



which are satisfied by almost all c and uj for all nonzero t. Proceeding as in the proof 
of Theorem I2.H 



|/| < CY,\M>t'\\^ - cos^,)-i(l - cos^,0"'(|l - e^(2c+(^+^'V)|-i + |i _ e' 



oo 

fcl,fc2 = l 
OO 

n'" \ ^ i.m-1 i,-2r i.-2r im+e im+e ( im+t . i.m+t\ 

^ 7 Aj-j^ 1^-^ i^<2 '^'2 V 1 2 / 

fcl ,^2 = 1 

OO oo 
= 2C"" ^ ^-{l+2{r-3m/2-£)) ^ ^-{l+2{r-m-e/2)) ^ ^ 



fci=l fc2=l 

provided we choose e > so small that r > 3m/2 + e. ■ 
6.2 Chaotic case 

Recall that Kc = 1 in definition (11.31) if and only if Mc{n) = an + o{n) where a > 0. 
Equivalently Dc{n) = an + o{n) with a > 0. We show that this is a sufficient condition 
for Kc = 1 via the correlation method. 

Theorem 6.2 Let c E (0, 2n). Suppose that Dc{n) = an + o{n) for some a > and 
that Kc is computed using the correlation method. Then = 1. 

Proof We claim that cov(,^, A) = -^an"^ + o(n^) and var(A) = j^a'^n'^ + o(n^). The 
result is then immediate. 

We verify the claim for cov(^,A). The verification for var(A) is similar. Write 

Dc{n) = an + e{n) where e{n) = o{n). Then 



(e, A) = lj2^iaj + e{j)) - (1 (1 Y.(aj + e^))) 

j=i j=i j=i 

1 1 " 1 1 " 

= -a(n + l)(2n + 1) + - " ^("^ + 1){ 2«^(^ + 1) + 

1 1 " 1 

-an^ + 0{n) + - je(j) - — (n + 1) e(j). 



12 ' ' n^-" 2n 
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Hence it remains to show that je(j) = o{n^) and — o{n^). 

Since e{n) = o{n), there is a constant C > such that |e(n)| < Cn for n > 1. 
Also, given e > 0, there exists uq > 1 such that \e{n)\ < en for all n > no. Choose ni 
such that Cn^/n\ < e. Then for all n > ni, 

^ n ^ no n ^ ^ n 

j=l j=l j=no+l ^ j=l 

< e + e = 2e, 

so that J2j=i^U) = o(n^). Similarly ^^^ije{j) = o{n^). ■ 



7 Discussion 

We have addressed the issue of validity of the 0-1 test as presented in |6l [7]. The 
original 0-1 test [5] included an equation for a phase variable 6{n + 1) = 6{n) + c + 
V o f"- rather than a constant "frequency" c. The extra equation for the phase driven 
by the observable made available theorems from ergodic theory on skew product 
systems |H [13], |15] . These theorems rely on the fact that for typical observables v the 
augmented system with the phase variable is mixing. 

In the modified version proposed in [6], the phase variable is not mixing, and the 
results from ergodic theory are not applicable anymore. Nevertheless, the modified 
test is more effective, particularly for systems with noise [B]. In this paper, we have 
verified that the modified test can be rigorously justified. Moreover, our theoretical 
results are stronger than the corresponding results mentioned in [5] for the original 
test. 

Our main results in this paper are that Kc = with probability one in the case of 
periodic or quasiperiodic dynamics, and that Kc = 1 with probability one for "suffi- 
ciently chaotic" dynamics. The latter includes dynamical systems with hyperbolicity, 
(including weakly mixing systems such as Pomeau-Manneville intermittency maps). 
In particular, nonuniform hyperbolicity assumptions combined with summable auto- 
correlations for the observable suffice to obtain Kc = 1. In the absence of hyperbol- 
icity, we still obtain Kc = 1 (with probability one) for observables with exponentially 
decaying autocorrelations. These results extend to systems that are mixing up to a 
cycle of finite length. 

We also made explicit the connection with power spectra: the test yields Kc = 1 
with probability one if and only if the power spectrum is well-defined and positive 
for almost all frequencies. The criteria above - exponential decay of autocorrelations 
or summable correlations plus hyperbolicity (up to a finite cycle) - are sufficient 
conditions for existence and positivity of the power spectrum. 
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There remains the question of whether typical smooth dynamical systems are 
either quasiperiodic or have power spectra that are defined and positive almost ev- 
erywhere. This is required for a complete justification of the test for chaos. Unfortu- 
nately the current understanding of dynamical systems is inadequate to answer this 
question, but all numerical studies so far indicate this to be the case. We leave it as a 
challenge to the skeptical reader to concoct a robust smooth example where the test 
fails! On the positive side, we showed in this paper that under a mild assumption on 
autocorrelations, slightly stronger than summable but much weaker than exponen- 
tial, we obtain either i^c = or Kc = 1 for each choice of c, though without invoking 
hyperbolicity we cannot rule out the possibility that both Kc = and Kc = I occur 
with positive probability. 

Our investigations of the validity of the test for chaos enabled us to construct an 
improved version of our test. The modification which amounts to using Dc{n) rather 
than the mean square displacement Mc{n) was shown to significantly improve the 
test in [7]. In addition, we showed in [7j that Kc is better computed by correlating 
the mean square displacement with linear growth rather than computing the log- 
log slope. In this paper, we have shown that our rigorous results apply also to the 
improved implementation of the test in [7]. 
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